Moving contact line of a volatile fluid 
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Interfacial flows close to a moving contact line are inherently multi-scale. The shape ol the 
interface and the flow at meso- and macroscopic scales inherit an apparent interface slope and 
a regularization length, both called after Voinov, from the dynamical processes at work at the 
microscopic level. Here, we solve this inner problem in the case of a volatile fluid at equilibrium with 
its vapor. The evaporative/condensation flux is then controlled by the dependence of the saturation 
temperature on interface curvature - the so-called Kelvin effect. We derive the dependencies of the 
Voinov angle and of the Voinov length as functions of the substrate temperature. The relevance of 
the predictions for experimental problems is finally discussed. 



The dynamics of a macroscopic solidplunging in a liq- 
uid bath [J 0] or withdrawn from it [j-HJ depends sen- 
sitively on its wetting properties i.e. on the intermolec- 
ular interactions at the nanoscopic scale. The motion of 
the contact line separating wet from dry regions is there- 
fore an inherently multi-scale problem. Amongst the im- 
portant consequences of the coupling between inner and 
outer scales (Fig. QJ, the speed at which a contact line 
can recede over a flat solid surface cannot exceed a crit- 
ical value, associated to a dynamical wetting transition 
which leads to the formation of a dewetting ridge 



of a V-shaped dewetting corner 
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or to the en- 
for detailed 



trainment of films 
reviews). In many applications, such as coating, imbi- 
bition of powders, immersion lithography or boiling-free 
heating, these entrainment phenomena are crucial limit- 
ing factors for industrial processes. 

Fig. Q] shows schematically the structure of the flow 
close to a moving contact line. Even for an infinitesi- 
mal velocity U, there exists a range of mesoscopic scales 
- roughly six decades - separating the microscopic scale 
from the macroscopic length L, in which the diverging 
viscous stress is balanced by a gradient of capillary pres- 
sure. This balance can be made quantitative in the lubri- 
cation approximation, for which the angles are assumed 
small, and which gives a third order differential equation 
for the interface profile h(x): 
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dx 3 



3r/U 



(1) 



where 77 is the liquid dynamic viscosity and 7 the surface 
tension; U is positive for an advancing contact line. This 
equation has an exact solution [l4| which reduces to the 
asymptotic form proposed by Voinov fl5| far from the 
contact line, but for x <C L: 



h'(x) 
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By is by definition the apparent contact angle in the 



FIG. 1. (Color online) Schematic showing a liquid vapor in- 
terface h(x) - here, a spreading drop (U > 0) on a cooled 
plate (AT < 0) - at different scales. The inner region, 
close to the moving contact line, is controlled by evapora- 
tion/condensation (top right). The slope changes from the 
Young angle By to the Voinov angle 8v across a scale lv given 
by the Kelvin length Ik- In a mesoscopic range of scales, 
the shape of the interface results from the balance between 
viscous friction and Laplace pressure gradient, resulting in a 
slope h'(x) which varies logarithmically in scale (Eq. [2]). 



static case (U — 0), which can be different from the 
Young angle 0y due to out of equilibrium processes tak- 
ing place at a microscopic scale. The Voinov length ly is 
also a quantity defined in the mesoscopic range of scales 
but inherited from the inner region, where the problem 
is regularized. The mesoscopic solution @ must also 
be matched at the macroscopic scale L to an outer so- 
lution where viscosity can usually be neglected. Fig. Q] 
features the case of a spreading drop or of a growing 
bubble but the outer matching problem has been solved 
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for manyother geometries including a gravity controlled 
bath [1,13, EH or a capillary ridge Q ■ 

Different models have been proposed to solve the mov- 
ing contact line paradox i.e. the singularity of Eq. (P) 
as x — > 0. The simplest regularization is obtained by 
imposing the Navier slip boundary condition based on 
a slip length £ s that can be expressed using a statisti- 
cal physics description of liquids [T3| and gas (HI. The 
Voinov angle is then Young's angle Qy and the Voinov 
length reads iy = 3£ s /(e 6y), where e is Euler's num- 
ber Jij. Alternative descriptions have been proposed, 
based on disjoining pressure (ly then scales on the Is- 
raelachvili's length (A/Gtt^/) 1 ^ 2 , where A is the Hamaker 
constant) and diffuse interface models [ly is then set by 
a diffusion length). Finally, when the substrate present 
heterogeneities, the contact line dynamics becomes in the 
inner layer a thermally activated process [l3|, EH . 

Here we consider the contact line motion of a volatile 
liquid in contact with an atmosphere of its pure vapor, 
2(3- HI and references therein. 



see 



We show that the 
lubrication equation is perfectly regular when evapora- 
tion/condensation processes are taken into account. The 
Voinov length £y and the Voinov angle By as functions of 
the substrate temperature constitute the central results 
of this Letter. The relevance of the theory for the case 
of a drop evaporating in air, usually assumed to be con- 
trolled by vapor diffusion [25l - l28f , will be discussed in the 
conclusion. 

Lubrication equations including Kelvin effect - Evap- 
oration/condensation process has first been proposed as 
a possible mechanism controlling the contact line motion 
at the molecular scale by Wayner [2l| and Pomeau [30j ]. 
In this Letter, we formalize this idea in a hydrodynami- 
cal framework, starting from the equation governing the 
evolution of the interface position h: 



dih + d x q = -j 



(3) 



where q is the hydrodynamic flow rate. The rate j at 
which a liquid evaporates is governed by the energy bal- 
ance at the liquid-gas interface. Assuming that the va- 
por pressure is fixed, the interfacial temperature T l de- 
pends on the interface curvature k ~ h"(x) according to 
Kelvin's law 
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where T s is the saturation temperature, C the latent heat 
and p the liquid density. In the lubrication approxima- 
tion, the temperature varies linearly across the liquid 
layer from the substrate temperature T s + AT, assumed 
to be imposed (Fig. [I]), to the interfacial temperature Ti. 
Neglecting the energy flux in the vapor phase, the evap- 
oration rate is controlled by the conductive energy flux 
across the liquid, 
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FIG. 2. Solution of the equations for a receding contact line 
(S = —0.02), with an overheating e = 0.1. The dotted line 
corresponds to the static case (5 = 0). The dashed line corre- 
sponds to the Voinov outer asymptotics ([2]). These solutions 
allow one to obtain the Voinov length ly and angle 6v- 

where k is the liquid heat conductivity. 

Starting from the Voinov law ([2]), we define the reduced 
capillary number [31(, using Qy as a characteristic slope: 



5 = 



ZnU 



(6) 



The characteristic length £k is obtained dimensionally 
by balancing the two fluxes driven by the interface cur- 
vature, namely the evaporation rate j and the divergence 
of the hydrodynamical flow rate q = jh 3 k' j (3n): 



£ K = 



\JZr}kT s 
~0 v ~p~£~ 



(7) 



We therefore make the solution dimcnsionless using 

h(x) = 6 v l K H (C) , C = x/£k, 

Under the lubrication approximation, the governing 
equations in the scaled variables then read 

H" = K, Kf = QH~ 3 - 5H~ 2 , Q' = (IC-ej/H, (8) 

where e is the superheating parameter, defined by 



y/3l]kT s AT 



(9) 



and Q is the dimensionless counterpart of (q + Uh). The 
fourth order differential equation (JSJ) must be comple- 
mented by the appropriate boundary conditions. We 
choose x = for the contact line position (so H(0) = 
0) and impose the slope H'(0) = Qy/Ov according to 
Young's law. As we look for regular solutions of the 
problem, the continuity of the temperature at the con- 
tact line requires /C(0) = e. Finally, to make the problem 
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compatible with the asymptotic expansion ([2|), one as- 
sumes a vanishing curvature far from the contact line: 
/C(oo) ->■ 0. 

Voinov angle - By definition, 9y is the interface slope 
for i > 1^ at vanishing capillary number S. In this limit, 
the outer boundary condition /C(oo) — >■ is equivalent to 
the constant slope condition H'(oo) = 1. The dotted line 
in Fig. [2] corresponds to a typical solution obtained nu- 
merically for (5 = Q. Overheating (e > 0) induces an 
evaporation flux that would diverge as at the con- 
tact line (cf. Eq. [SJ), if not balanced by Kelvin's effect. 
The induced liquid flow towards the contact line is ac- 
companied by a capillary pressure gradient: the resulting 
interface curvature leads to an apparent angle 9y larger 
than 9y . The cross-over from 9y to 6y takes place at 
the microscopic scale, for £ of the order unity. The ratio 
9y/9y is reported in Fig.[3k as a function of e. Following 
[23(, we perform a linear expansion of the solution in e 
for (5 = 0, writing H = Hq = C,+eH e + 0(e 2 ). Linearizing 
Eq. [8j one obtains a differential equation on K, e = H' e ', 
which reads (C 3 /C^)' — /C e /£ = — 1/C The solution verify- 
ing the boundary conditions involves the modified Bcsscl 
function of the first order K\\ K. e = 1 — C,~ 1 K\ (C 1 )- 
Integrating /C e from oo to 0, one obtains H' e (0) = —tt/2, 
which gives the expansion for the Voinov angle 



? Y /8 v = l-(ir/2)e + 0(e 2 ), 



(10) 



shown in dashed line in Fig. [3k. 

The most important feature of the curve 9y/9y(e) is 
the existence of a critical value e c ~ 0.297 of the over- 
heating parameter - note that the linear approximation 
(fTOf overestimates e c by a factor ~ 2. In the limit e — > e c , 
9v becomes much larger than 9y . Then, Kelvin effect is 
just balanced by the maximal available capillary force, so 
that H' goes from to 1 [34[. The equation e = e c gives 
the large AT asymptotic expression of the Voinov angle 
9y, which does not depend any longer on 9y: 



{ V3ijkT s AT 
V £ c7 Ts 



1/3 



(11) 



One may expect this asymptotic regime to be relevant 
close to the gas-liquid critical point, in particular to the 
description of boiling [32. 33 1. 

Voinov length - We now consider a contact line mov- 
ing at a velocity U. We linearize the governing equations 
with respect to (5, around the solution Hq obtained for 
5 = [34| . A typical solution is shown in Fig. [2] (solid 
line) together with the asymptotic expansions around 
C — > (static solution obtained for (5 = 0, dotted line) and 
C — > co (Voinov expansion, dashed line). It shows that a 
perfectly regular solution is obtained, in spite of the no- 
slip boundary condition imposed at the solid-liquid inter- 
face: at a scale smaller than £k, the interface advances by 
the curvature driven condensation (or recedes by evapo- 
ration). How can a contact line advance (even at e = 0) 
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FIG. 3. (a) Ratio of Young angle to Voinov angle as a function 
of overheating parameter e, determined numerically (solid 
line). The dashed line is the analytical expansion (IIOI) . 



in the absence of any regularizing mechanism leading to 
a slip of the fluid at the boundary? Consider a perfect 
wedge initially at rest at the Young angle. Imposing an 
hydrodynamics flux towards the contact line leads to an 
increase of the apparent contact angle. However, 9y re- 
mains to be the true contact angle at the molecular scale 
so that a positive curvature k appears at a small scale, 
which induces a condensation by Kelvin effect. As a con- 
sequence, the contact line advances although the liquid 
velocity vanishes at the contact line: the phase transition 
flux j(0) = —U9y balances exactly that induced by the 
contact line motion. 

At distances much larger than the scale £y, one re- 
covers as expected the Voinov solution H'(() ~ 1 + 
S log(C^K /£v ) + O ((5 2 ) . The Voinov length £ v is obtained 
from the matching to this outer expansion, as shown ge- 
ometrically in Fig. [5] (intersection between the dashed 
line and the horizontal line H' = 1). Fig. [3] shows the 
dependence of £y on the overheating parameter e. As 
expected from dimensional analysis, £y is on the order 
of the Kelvin length £k ■ The ratio £y / £ k turns out to 
increase with the overheating parameter, from ~ 1.32 at 
e = to ~ 3.00 at e = e c . 
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Concluding remarks - As long as one aims to address 
a macroscopic problem involving a moving contact line, 
the only quantities inherited from the inner molecular- 
scale region are the Voinov length and the Voinov an- 
gle. For a volatile fluid with a vanishing slip length, we 
have shown here that the Voinov length is set by the 
Kelvin length. However, a true fluid presents both slip at 
the solid/liquid interface and evaporation/condensation 
at the solid/vapor interface. The theory developed in 
this Letter is applicable if the Voinov length produced 
by Kelvin effect is larger than that induced by the slip 
length i.e. if the condition dy Ik/^s > 1 is fulfilled. The 
product 6 V Ik depends only on the liquid properties and 
ranges from 0.3 nm for water and methanol to 1 nm for 
alkanes and refrigerants such as ammonia or fluorocar- 
bon. It can be even larger, for fluids like glycerol or sili- 
con oils whose large viscosities are due to glassy effects. 
As the slip length £ s is around two molecular sizes when 
Qy < 7r/2, a good proxy of 9y £kI^s is v l : evaporation 
is indeed the regularizing mechanism for any contact line 
problem, in the limit of low contact angles. 

The theory developed here is directly applicable to the 
microscopic description of boiling, as nucleating bubbles 
are constituted of pure vapor. It may help resolving the 
demanding problem of the boiling crisis 33j. However, a 
second evaporative effect has been neglected here: when 
a fluid molecule joins the vapor phase, its momentum is 
lost by the liquid. A normal stress p 2 j 2 /p v is therefore 
exerted on the interface (p v <C p is the vapor density). 
Balancing this effect with capillary pressure 7 /Ik, one 
deduces that the vapor recoil effect is dominant when 
the contact line velocity U is much larger than the recoil 
velocity U r = (-fC)^ 2 (p v / p) 1 / 2 {SnkT,)- 1 ^. U r is very 
large (~ 10 m/s for a fluorocarbon) so that vapor recoil 
pressure can usually be neglected. 

Is Kelvin's effect relevant to the coffee stain prob- 
lem [i^] where a drop evaporates in air? The rate of 
evaporation j in such a case must satisfy simultaneously 
three conditions: (i) the kinetic equation based on Hertz- 
Knudsen law, (ii) the conservation of mass, controlled by 
the diffusion of vapor and (iii) the energy conservation 
at the interface. The problem reduces to that consid- 
ered here if the first two processes are fast enough, i.e. if 
the evaporation rate predicted by Eq. is the limiting 
factor, (i) The kinetics of evaporation involves thermal 
velocities of a fraction of the speed of sound in air, much 
larger than observed evaporation rates, (ii) The evap- 
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oration rate j predicted by the diffusion alone 
diverges at the contact line as a; -1 / 2 , (iii) In conclu- 
sion, the evaporation rate at the microscopic scale must 
be generically controlled by the energy brought by ther- 
mal conduction to the liquid/air interface. Kelvin's effect 
shall therefore regularize both the stress and the evapo- 
rative flux singularities induced by the wedge geometry 
of a contact line 351. 
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